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■ We establish the critical line of the one-component $ 4 (or Landau-Ginzburg) model on the 

simple cubic lattice in three dimensions. Our study is performed in the framework of the non- 
perturbative renormalization group in the local potential approximation. Soft as well as ultra- 
sharp infra-red regulators are both considered. While the latter gives poor results, the critical 
line given by the soft cut-off compares well with the Monte Carlo simulations data of Hasenbusch 
(J. Phys. A : Math. Gen. 32 (1999) 4851) with a relative error of, at worst, ~ 3.10" 3 on published 
points (critical parameters) of this line. 
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I. INTRODUCTION 



The non perturbative renormalization group (NPRG) approach initiated by Wetterich et al. 1-3] has proved its 
ability to describe both universal and non universal quantities for various models of statistical and condensed matter 
!_? 1 physics near criticality. It has been extended and adapted to lattice models recently [4[ and applied successfully to the 
three-dimensional (3D) Ising, XY, and Heisenberg model Here we apply the lattice NPRG to the one-component 
<I> 4 model in 3D on a simple cubic lattice. 

We work in the framework of the local potential approximation (LPA) and consider both a sharp and a smooth 
infra-red cut-off (or regulator). As in Refs. [a, 0] the flow equations are integrated out for the so-called threshold 
functions rather than for the potential. The resulting flow equations turn out to be quasi-linear parabolic partial 
differential equations (PDE) for which several stable and unconditionally convergent numerical algorithms have been 
developed by mathematicians § . We made use of the algorithm of Douglas- Jones [j| @ to solve the NPRG flow 
. equations both above and below the critical temperature; this yields an easy and precise determination of the critical 
point. The critical line of the model is obtained for a large range of parameters and compared with the Monte Carlo 
t^J- I (MC) data of Ref. [Icj. 

t-H i Our paper is organized as follows : In Sec. [Hlwegive a short review of lattice field theory and of recent advances 
in the application of the version of Wetterich [l], [2} of the NPRG to lattice systems [3, @- Then, in Sec. IIII1 we 
detail the LPA approach which constitutes the simplest non-perturbative approximation to solve the flow equations. 
Two convenient regulators are introduced and the corresponding flow equations are derived. Mathematically, the 
flow equations for the potential are non-linear parabolic PDE for which the correct initial conditions need a thorough 
discussion. A numerical solution of these equations requires a change of variables explained in Sec. IIVI After this 
transformation, the emergent equations turn out to be quasi-linear parabolic PDE and can thus be solved easily on a 
computer making use of powerful algorithms. Numerical results for the critical line of the one-component $ 4 model 
on the simp le cubic lattice in three dimensions are reported in Sec. IVland compared with the Monte Carlo (MC) data 
of Ref. [10| • The sharp cut-off regulator leads to poor results while the soft cut-off one reproduces the MC data with 



1 a relative precision of a few 10 . We conclude in Sec. [VTJ 

II. PROLEGOMENA 
A. Lattice $ 4 model 

We consider a lattice field theory defined on a D-dimensional (hyper)cubic lattice. The action is given by [TTI j 
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where a is the lattice constant, and {r} denotes the N sites of the lattice. For simplicity we limit ourselves to a 
one-component real field <p r and a simple cubic (SC) lattice, </? q = a D X){ r } e ~ lrt Vr is the Fourier transform of the 
field and the ./V momentum {q} are restricted to the first Brillouin zone ] — ir/a,ir /a] D . In the thermodynamic limit 
(a fixed, N — > oo), 
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The local potential Uq is defined such that eo(q = 0) = and henceforth we adopt the Landau- Ginzburg polynomial 
form Uo(cp) — (r/2) ip 2 + (g/4!) ip 4 . The spectrum e (q) accounts for next-neighbor interactions; for a SC lattice it 
reads as 



D 



eo(q) = (l-cos(^a)) , 

a n=l 



(3) 



where the dimensionless parameter £ rules the amplitude of the spectrum, eo (q) ~ £q 2 for q — > and max q eo (q) = 
£ max _ (4£)£)/ a 2 . Note that the dimension of the field is [ip] = D/2 — 1 and that, in the thermodynamic limit, the 
physics of the model depends only upon the two dimensionless parameters r = ra 2 /£ and g = ga( 4 ~ Z3 - ) /£ 2 



Another way of writing the action ([!]), which is useful for numerical investigations, i 



and g = 

,is M 



{n} L 
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where the D unit vectors e M constitute an orthogonal basis set for MP . The field ip and the parameters (k, A) are all 
dimensionless awnd are related to the bare field ip and dimensionless parameters (r, g) through the relations 

"0n = \j a / 2 ^ 1 ip r with r = an , (5a) 



1 - 2A 
r = 2D 

K 
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(5b) 
(5c) 



To close this section let us recall some elementary definitions and results concerning lattice field theory [Tij . The 
physics (thermodynamic and correlation functions) of the model is coded in the partition function 



Z[h] 



Vtpexp(-S[(p] + {h\<p)) , 



where the dimensionless measure is given by 



np = H [< 



D/2-1 



dipr 



h is an external lattice field and the scalar product in ([6]) is defined as 

(h\ip) = a D h r (p T ■ 



(6) 



(7) 



(8) 



The order parameter is given by 



(9) 
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where the Helmholtz free energy W [h] =lnZ [h] which is a convex function of the N variables {h r } is the generator 
of the connected correlation functions G^ n \vi . . .r n ) = a~ nD (d j dh T1 ) . . . (d / dh rri )W [h]. The Legendre transform of 
W [h], i. e. the Gibbs free energy, will be denoted unusually by T [(f)] and reads as 

T[4>] = {h\<t>)-W[h] . (10) 

— —M 

r [cp\ is a convex function of the N field variables {4> T } and the generator of the vertex functions T (r x . . . r„) = 

a~ nD (d/dcj) ri ) . . . {d/d(f> rn )T [(f)]. Finally, as well known, the matrix r 2 (ri, is the inverse of matrix (ri, r^). 



B. Lattice NPRG 

To implement the lattice NPRGprocedure we follow the suggestion of Dupuis et al. in Ref. 0, Q which extends to 
the lattice the ideas of Wetterich for the continuous version (a — > 0) of the model. We thus add to the action (fTJ) 
the regulator term 

A5fcM = ^A^^- qiife(q) ^ q - (U) 
{q} 

where Rk (q) is positive-definite, has the dimension [Rk] = 2 and acts as a q dependent mass term. The regulator 
Rk (q) is chosen in such a way that it acts as an infra-red (IR) cut-off which leaves the high-momentum modes 
unaffected and gives a mass to the low-energy ones. Roughly Rk (q) ~ for ||q|| > k and Rk (q) ~ Zkk 2 for ||q|| < k. 
The scale k in momentum space varies from A ~ a" 1 , some undefined microscopic scale of the model yet to be 
precised, to k — the macroscopic scale. To each scale k corresponds a fc-system defined by its microscopic action 
Sk [f] = S [<p] + A5fc [ip] . We denote its partition function by Zk [h] , its Gibbs free energy by Tk [4>] , etc. 

For technical reasons that should become clear below, we are rather interested in the so-called average effective 
action Tk [(f)] which was introduced by Wetterich [l| and is defined as a modified Legendre transform of Wk [h] which 
includes the explicit substraction of ASk [4>] 0> *• e - 

T k [(f)] = Tk [<f>] - ASk [<f>] ■ (12) 

Note that the functional Tk [(f)] is not necessarily convex by contrast with T [(f)] which is the true Gibbs free energy of 
the "k"- system. It satisfies the exact flow equation [l|- 



1 



dk Tk [<i>] = - J2 d * R * (q) r i 2) + R k _ .. ■ ' I - • 



q 



( 2) 

Note that Eq. (fl3|) is an extremely complicated equation since the vertex function 1^ (q, — q), which is the Fourier 
transform of the second-order functional derivative of T [(f)] with respect to the classical field 0, depends functionally 
upon (f>. For an homogeneous configuration of the field (p T = (f> we have, on the one hand, Tk [4>] — Na D Uk{4>) where 
the potential Uk{4>) is a simple function of the field <f> and, on the orther hand, the conservation of momentum at each 
vertex which implies, with the usual abusive notation, r^(q, — q) = Na D T k 2 \q) from which follows : 



a tt fA\ 11 dkRk (l) ha \ 

^)^^E r ( 2)(q)+i?fc(qr ("a) 

i r d k R k ( q ) 



2 7 q ri 2) (q)+i?,(q) ' 

where the second line (|14bl) is valid in the thermodynamic limit (a fixed, N — > oo). 
We can give a formal solution of fp~3|) as [H, Q 

e*p(-r fc [fl) = Jv<f>exp(-S[<p] + (<p-<l>\^^)-^-4>\Rk\<p-'l>)j (15) 

which gives Tk [(f)] implicitly. Eq. (|T5|) allows us to precise the initial conditions. The initial value k = A of the 
momentum scale k is chosen such that i?A(q) ~ oo for all values of q hence, since exp(— 1/2 (xI-^aIx)) ^ £[x])> where 
S[x] is the Dirac functional, it follows from (TTSl) that T\[4>] — S[<f>], Physically it means that all fluctuations are frozen 
and the mean- field theory becomes exact. When the running momentum goes from k = A to k = all the modes tpq 
are integrated out progressively and the effective average action evolves from its microscopic limit = <S[</>] to its 

final macroscopic expression Tk=o[(f)] — r [(/>]. The choice of initial momentum A which depends on the choice of the 
regulator Rk will be made more explicit in Sec. Mil 
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III. LOCAL POTENTIAL APPROXIMATION 



An increasingly popular way to solve the flow eq. (| L3[) is to make an ansatz on the functional form of the effective 
average action r^[0]. In the lattice LPA one neglects the renormalization of the spectrum and assume the local 
form 

(LPA ansatz) T k [4>] = — L ^ ^_ q£o ( q ) q + a D £ U k (cf> r ) . (16) 

{q} M 

For a uniform configuration of the classical field <fi r = cf> and, in the thermodynamic limit, the flow equation ()14bp 
becomes : 

d k U k ((j>) = - / - , 17 

2 J q eo(q) + Rk (q) + t/ fc w 

where U k ((f)) denotes the second-order derivation of Uk((f>) with respect to the order parameter <j>. Eq. ((TT)) is a 
non-linear parabolic PDE. It must be supplemented by an initial condition (see Sec. MI Al and (|III B^ ). Moreover, for 

a numerical resolution of ()17|) . boundary conditions for the potential Uk or one of its derivatives uj^\ must also be 
specified for some maximum value of the field ±^> max (see Sect. IIV[) . Initial and boundary conditions depend on the 
choice of the regulator Rk and, in this paper, we will consider two possibilities for Rk- 

A. Litim-Machado-Dupuis (LMD) regulator 

In Ref. |5| Machado and Dupuis consider 

R k (q) = [e fc - e (q)] 6 [e fe - e (q)] , (18) 

which is adapted from Ref. [13] to the lattice case. In Eq. (TT5)) . = £ fc 2 and 0(a;) is the step function. This regulator 
i?fc(q) leaves the high-momentum modes (eo(q) > £fc) unaffected and ascribes a mass to the low-energy ones. The 
effective spectrum of the k-system is obviously e£ ff (q) = eo(q) + i?fc(q). We note that for k > fc max , where fc max is 
defined by e/c max = to 13 *, %■ e. fc max = y/4D/a 2 , the effective spectrum e| ff (q) = €k does not depends on q and we 
deal with a theory of N independent sites. 

With that choice the flow equation (fl7|) takes then a simple synthetic form 

d t Uk = ^-77 M{e k ) (19) 

where the RG time "i" is defined by k = Ae - *, so that dt = —kdk, and 

N(e k ) = f 0(e fc - eo(q)) • (20) 
denotes the (normalized) number of states. It will prove convenient to introduce also the density of states 

V(e k ) = I 6(e k - e (q)) , (21) 

so that 

Af(e) = I de V(e) . (22) 



q 



These two functions T>(E) and N(E) are obviously related to the lattice Green function which, for a SC lattice, 
reads QJ, 03 

= i f dgx • • . f d to . (23 ) 

Jo Jo t-}J^ = i cos(g M ) 
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FIG. 1: Density of states T>{E) and number of states M{E) of the simple 3D cubic lattice 



With the remark that for r] — > 0, l/(r + irf) = V(1/t) + i-k5(t), the comparison of eqs (|2"Tj) and (|2"3"f reveals that 



25(e) 



-3G(r) , 



(24) 



2£ 7T 

with r = D — (a 2 /2)(e/£). Note that the interval of the spectrum < e& < e™ ax corresponds to the interval 
— Z? < r < D for the auxiliary variable r. In the case of D = 3 the imaginary part of the Green function G(r) is 
given by fl4| 



1< r < 3 



< r < 1 



3GM = ^1 



1/2 > 



x 2 - 1 



(25) 



(26) 



where x = 2/(r — cos(x)) and K (y) is the complete elliptic integral of the second kind 

K{y)= I*'' deil-y^^ie))- 1 ' 2 . 
Jo 

One can point out the properties 3G(r) = 3G(— r) for —3 < r < 3 which implies that T>(e) = 2?(e™ ax — e) and 
j^max _ e ) = i _ 7V( e ) for < e < eJJ lax . 

We made use of the relations ([23)1 to evaluate numerically T>(E) a.nd J\f(E), their graphs are displayed in Fig ([1]). The 
tiny wiggles in the central part of T>(e) are actual and could not be avoided, they reveal the difficulty to compute this 
function with the highest numerical precision. A numerical filter can be used to suppress the numerical fluctuations 
in T>(E) and N(E) before the latter is injected in the flow equations (Q2 



1. The local regime : A > k > fc ma x 

Clearly one can distinguish 2 different regimes in the flow. For larger than the gap of the spectrum, tu = £fc 2 > 
e™ ax , i. e. k > fc max = y/AD/a 2 , we have already pointed out that the "effective" spectrum of the k-system e| ff - = e\~ 
does not depends on q. We deal with a local theory for which the partition function Zk [h] — Y[ r Zk{h r ) is a product 
of one-site partition functions with 

Zk(h) = J dip exp (-U (cp) - \fW 2 + h<pj . (27) 



-2-1012 
3> 



FIG. 2: Local potential of the "I? 4 model on a 3D simple cubic lattice at g = 1, and r = —0.13. Dotted line : mean field 
approximation (A = 5000a -1 ), Dashed line : local theory at k — fc ma x, Solid line : Uk=o in the LPA approximation with LMD 
regulator. 



It follows from this remark that the effective average action has also a local form 

Tk M = ^ ^ (q) 0q + a ° £ u * c ' ^ > (28) 

{q} 

where the effective average potential U k oc ' is given implicitly by 

exp (-Ul oc - (4>)) = J + °° dip exp (-UoM + U£°- ' (<t>)(<p - <P) - \t k (ip - 4>?^j , (29) 

as follows from eqs (|15p and (|2"T|) . 

Two remarks are now in order. Firstly the choice A = oo implies [7a = Uq since we can replace the gaussian 
exp(— l/2e\((p — 4>) 2 ) by a delta function S(ip — <f>) in eq. (f29|) . Secondly in the range A > k > /c max the potential Uk 
satisfies the exact flow equation 

d t U k = ^-r , (30) 

as shown in the appendix. Note that for A > k > fc max we have T>(ek) = 1 and thus the LPA flow equation p^|) 
is exact for a local theory. This point has been checked numerically in Ref.Q and is proved mathematically in the 
appendix. The initial condition of the flow can thus be chosen 

• either A = oo and Ua = Uq (Mean field theory) 



• or A = fc max = y^AD/a? and U A = U k oc ; 



In the second case the local partition function (|27l) must be computed numerically. It turns out that integrating out 
the flow equation (|30j) with the mean field initial condition (numerically with a large value of A, e. g. A ~ 5000a -1 ) 
gives more accurate results than the direct calculation and manipulations of z(h) which involve too large arguments 
in the exponentials. 

We exemplify this discussion in Fig. [2] where we display, for the $ 4 model at r = —0.13 and g = 1 (i. e. in the 
ordered phase), U\((f>) — Uq(4>) (mean field approximation), J7fe max (<^) (local theory) and the renormalizcd potential 
Uk(4>) at k = obtained after integration of the flow equation (jT9")l (the figure also illustrates the passage to convexity 
: note the flat part of Uk=o{(t>) m the range (—4>o, <fro))- 
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2. The non-trivial regime : fc max > k > 



In this range of k the flow is non-trivial and must be integrated out numerically. For convenience we rewrite f| 19[) 

as 

8 t U k = -N{e k )C LMD {u k ) , (31) 
where uj k ((j)) = U k ((f)) /e k is a dimensionless renormalized susceptibility and 

C LMD (x) = T ^, (32) 

is the threshold function [2j which takes a very simple expression with the LMD regulator. 
Note that in the limit k — > 

q 



Af(e) ~ / 6(fc 2 - q 2 ) 
i 

1 



dq q D - X S, 



D 



{^) D Jo 

k D , (33) 



D 

where Sr> = 2n D / 2 /T(D/2) is the surface of the D-dimcnsional hypersphere and v^ 1 = 2 D+1 n D ^ 2 T(D/2) a numerical 
factor. Therefore in this limit the flow eq. p^|) reduces to 



dt U k = k D C LMD {u k ) , (34) 

which is identical to the LPA flow equation with Litim's regulator for the continuous (off-lattice) theory 0,0, EH- We 
conclude that the lattice (cf. eq. (flU)) ) and off-lattice (cf. eq. (JM))) how equations have the same asymptotic properties 
for k — > 0. The adimensionned versions of these equations share thus the same fixed points and the same critical 
exponents. All these quantities have been computed with the highest numerical accuracy in D — 3 dimensions, see 
e.g. Refs [7J, lla, LL6J . Recall that in the LPA, Fisher's exponent rj = in all dimensions of space (no renormalization of 
the spectrum) and that the other critical exponents are non-trivial and differ from the exact ones by a few per cents 
in D = 3. 



B. Ultra-sharp regulator 

The ultra sharp cut-off (USCO) was first introduced by Wegner-Houghton [l7| and considered by many authors in 
different NPRG studies of the continuous (off-lattice) $ 4 model |18l-[20|. In its lattice version, it also yields simple 
flow-equations. In this case, adapting the definition of Wetterich the regulator is defined as 

ifc(q) = Ze fc 9 [e fc - e (q)] , (35) 

where the constant Z is ultimately set to +00 [2|. In order to deal with the discontinuity of R k (q), we first introduce 
a smoothened version Q e (x) of the step function which varies mildly from to 1 in the interval (— e/2, +e/2). Let 
5 e (x) = d x Q e (x) be the smoothened version of the Dirac generalized function. Then the flow equation (fT7| takes the 
form : 



8 U (6) = - ( Ze kd k e k 6 e (e k - e (q)) If Zd k e k 9 £ (e fc - e (q)) 

k kW 2y q eo(q) + Ze fc e e (e fc -eo(q)) + C^'(0) + 2 7 q e (q) + Ze fe + [/£(«/>) 1 ' 

The limit e — > for the ill-defined first term in the r.h.s. of (|36]) can be taken by making use of an extension of a 
lemma due Morris felil which states that, for e — > 



limS e (e k - eo(q))/(9 6 (e fe - e (q),q)) = <5(e fe - e (q)) / dt f(t,q) 



(37) 
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provided that the function f(Q e (q, k), k) is continuous at k = q in the limit e — > 0, which is the case here. This yields 



dkU k {4>) = ^(d k e k )V(e k )\n 



f-k + U k ((f)) + Ze k 



+ ^(d k \ne k )Af(e k ) + 0(Z- 1 ) . (38) 



The last step is to take the limit Z — > oo of the above eq. with the final result 

d t U k = -e k V{e k )C usco {u> k ) , (39) 

where u) k (cf>) = U k (4>) / e k and the USCO threshold function reads 

C usco (x) = -ln(l + .T) . (40) 

Note that U k (4>) is defined up to an additive constant, i. e. independent of the field <fi, which was discarded from 
Eq. dSHD- 

Since, in the limit k — > 

2>(e)~ / <5(fc 2 -q 2 )/e 



~ 2 -fk°-* , (41) 

then, the flow eq. (|39"1) reduces to 

d t U k = -2v D k D C usco {u k ) (fc^O), (42) 

which is the LPA flow equation with an USCO regulator for the continuous (off-lattice) theory 0, [l8l - l20| . The 
fixed-points and critical exponents of the lattice and off-lattice versions of the 2 theories are thus identical. 

We now discuss the problem of the initial conditions. The USCO regulator R k (q) = Z = oo for all k > fc max = 
\J ADja 2 . It transpires from the discussion of sect. (|T8|) that the mean-field solution U k (4>) — Uq(4>) should be solution 
of the flow-equation (|39| for al A > k > fc max . Indeed for k > fc max we obviously have V(e k ) = from which d t U k = 
follows. Any A > fc max (with — oo) can be kept as a valid initial condition since the MF solution U\ — Uo does 
not evolves in the range A > k > A: max . 

A last remark is in order. The initial condition £4 max (0) = ( r /2)<^ 2 + (.g/4!)</> 4 yields a solution of the flow 
equation (|39|) for k < fc max only if 1 + Wfe max (4>) > for all <fi. Therefore the LPA with ultra sharp cut-off is defined 
for a negative r only if |r| < fc^ ax = AD /a 2 ; otherwise the flow equation has no solution. 

IV. A CRUCIAL CHANGE OF VARIABLES 

We pointed out in section IIII Al and IIII Bl that in the asymptotic limit k — > the lattice and off-lattice LPA 
flow equations bear the same form, both with LMD and USCO regulators. In the ordered phase this behavior is 
singular and has been studied at length in Refs @,0|- Briefly, in the limit k — > 0, uj k (4>) = U' k ' (4>) / (^k 2 ) —> —1 for 
— 4>o{k) < 4> < (f>o(k) where 4>o(k) is a precursor of the spontaneous magnetization 4>o — hm^o (f>o(k). It follows 
that the threshold functions C diverge in this interval as k 2 ~ D (for D > 2). This yields a universal behavior 
£(cf>)/£(4> = 0) = 1 — <j> 2 /cf)Q. Moreover, as a consequence, U k ((f>) becomes convex as k — > 0, in particular it becomes 
constant for — (f> < <j) < +(j} (see e. g. Fig. 

The divergence of the threshold functions makes impossible to obtain numerical solution of the non-linear PDE (fT9|) 
and (|39[) in the ordered phase, we really deal with stiff equations In order to remove stiffness, one is led to make the 
change of variables U k (M) L k (M) — C[uj k (M) = U k (M)/e k ]. We then obtain the equations 



(LMD) L'i{cj>) 



M{e k ) 



1 

1 



L k {4>) 



' dtL k (4>), (43a) 



M{e k ) L k (<f>Y 



(USCO) L k \<t>) = -f- [exp(-L fc fo)) - 1] + — cxp(-L fc (0)) d t L k {ct>) , (43b) 

where k = Ae _ *. 

By contradistinction with eqs (l3"Tj) and (|3^|) these quasi-linear parabolic PDE can easily be integrated out. As in 
Refs. [11,0] we made use of the fully implicit predictor-corrector algorithm of Douglas- Jones 0]. This algorithm is 
unconditionally stable and convergent and introduces an error of 0((Ai) 2 ) + 0((A(j>) 2 ) (At and A</> discrete RG time 
and field steps respectively) and can be used below and above the critical point as well. 

In order to solve eqs. (|43[) numerically one must precise the initial and boundary conditions. 
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t 

FIG. 3: U k (M = 0)/ejb as a function of the RG time "t" in the LPA approximation with LMD regulator. A = 1.145 is fixed 
and k c is obtained by dichotomy on k. Dashed lines : k > « c - Solid lines : k < k c - 

• (i) Initial conditions : they were discussed thoroughly in Sec. IIIII we have just to transpose this discussion to 
the threshold functions. For the LMD regulator one has, at t = 0, L kmax (<t>) = C LMD (U l h °^ {(f)) /e^ x ) for all 
\4>\ < '/'max, where m ax is the largest value of the field. In practice the local approximation Ui oc ((j))) and its 
derivatives with respect to the field are computed by integrating the exact flow ([29]) from its MF expression at 
some large A. For the USCO regulator one has, at t = 0, £ fcmax (» = C LMD {U'^ 0))/e™ ax ) for all |0| < max , i. 
e. one retains the mean field approximation of the potential. 

• (ii) Boundary conditions : for the LMD regulator we adopted Lfe(±^ max ) = £ LAID (Uj? c (0 m ax)/efe) for all k 
and for the USCO regulator Lfc(±</> max ) = C LMD (Uq (</> ma x)/efc) for all k 0j- It amounts to keep the first term 
in the hopping parameter (k) or loop expansions of Tk(4>) respectively, which is a reasonnable assumption at 
large fields. 

In order to determine the critical point one proceeds by dichotomy : for instance g is fixed and one varies r. The 
renormalized coupling constant uq = U k (M = 0)/ek in the limit k — > discriminates the state : for r > r c (g), 
uq — > oo, and for r < r c (g), uq — > —I. Alternatively one can fix A and vary k; an example is given in fig. El 

V. NUMERICAL RESULTS 

We solved Eqs (|43|) with the Douglas- Jones algorithm [9] in D = 3 dimensions of space. To fix the ideas we used 
for most our numerical experiments At = 10 , a maximum of N t = 1.810 5 time steps, A0 = 2.1CT 4 and = 15000 
field steps (i. e. max = 3.). For the LMD regulator the initial momentum was A = 5000a -1 . We checked that these 
values of the parameters give at least 7 stable figures for r c (g). Most our numerical studies were made by fixing the 
value of parameter g and varying r in order to determine its critical value f c {g) by dichotomy. 

Our data for the critical line f c (g) are given in table Q] for the USCO regulator, the eqs. (|43|) has no solution for 
r c < —12 which is a severe drawback. The data for the LMD smooth cut-off are displayed in table [Til i n this case 
the LPA equations admit solutions for all values of g and we stopped arbitrary our investigations at g = 1000. All 
these data are also displayed with the variables (A, k) in Fig. @] in order to be compared with the MC data that 
Hasenbusch obtained for several points [lol ]. As apparent in Fig. [4] the theoretical predictions of the LPA with an 
USCO regulator (crosses) are poor as soon as A > 0.5. We interpret this failure as a consequence of the use of a 
mean-field initial condition at A = fc max which turns out to be a bad approximation of the local theory at high values 
of A. By contrast a very good overall agreement between the Monte Carlo (MC) Data and the predictions of the LPA 
with LMD regulator is observed. This confirm the conclusions of Machado and Dupuis in Ref. [5j who obtained also 
such a good agreement in the case of the 3D Ising, XY and Heisenberg models. 

A more stringent comparison is made in table Mil where, for all the A considered in Ref. [Tcj, the critical k c (X) was 
obtained by dichotomy on k. The maximum relative error of the LPA-LMD theory can be seen not to exceed ~ 3.10 -3 
for the considered range of parameters. In the case of Ising, XY and Heisenberg models Machado and Dupuis report 
errors which are significantly higher, i. e. of the order of a few per cents, than the ones we obtained for the $ 4 





FIG. 4: Critical line k c (A) of the <E> 4 model on a 3D simple cubic lattice. Open circles : MC data of Ref. [13]; crosses : LPA 
with USCO regulator; solid circles : LPA with LMD regulator, a dashed line joints the points as a guide-line for the eyes. 
Uncertainties are smaller than the sizes of the symbols. 

model. However these authors used a standard explicit Euler integration scheme for the non-linear PDE for Uk which 
yields to stop the flow before its scaling limit k — > in the ordered phase We suggest that solving instead the 
quasi-linear parabolic equations satisfied by the threshold functions Lk could perhaps change the evaluation of the 
critical parameters, yielding a still better status for the LPA predictions. This point should be checked. 



VI. CONCLUSION 



In this paper we have computed the critical line of the $ 4 one-component model on the simple cubic lattice in three 
dimensions in the framework of the NPRG within the LPA approximation. We have considered both a sharp and a 
smooth regulator. The flow equations have been solved for the threshold functions rather than for the potential. This 
trick allows to obtain numerical solutions in the ordered phase where the PDE for the potential are stiff and fail to 
converge. 

A dichotomy process based on the generically different asymptotic behaviors of the adimensionned susceptibility 
U k (</) = 0)/fc 2 in zero field, below and above the critical point, yields a very precise determination of the (non- universal) 
critical parameters. 

The LPA with a sharp cut-off regulator must be supplemented with mean-field initial conditions at A = fc max which 
restricts the solution of the equation to a small domain of the (r, g) plane. Moreover, even in this restricted domain, 
the critical parameters are in poor agreement with the MC data of Ref. [l(|. When the smooth LMD regulator is 
considered, the PDE must be supplemented either with mean-field initial conditions at A = oo or by the exact local 
expression of the potential at any A > fc max . The various possible initial local conditions obey an exact flow equation 
which coincides with the LPA-LMD theory. The LPA with LMD cut-off gives surprisingly good estimates of the 
critical parameters of the lattice $ 4 model, the maximum deviation with MC data being of ~ 3.10~ 3 for the states 
which we considered. 

Extension of the present study to several other lattices seem feasible as well as the extension to vectorial O(N) 
models. 
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Appendix: Exact flow equation for the local theory 

We noted in the text that in the range A > k > fc max the average effective action is local for the LMD regulator. 
The local potential U k (4>) satisfies exactly to Eq. (j2"5)) and we rewrite this relation as 

exp (-U k (<f>) + U' k (cf))<t> + \tk4> 2 ^J = J dip exp (-U (ip) - ^e kV 2 + (u' k ((f>) + e k <f>)ip 

= z k (h) , (A.l) 

where z k (h) is the one-site partition function (|27[) and h = U k (4>) + e k (f) an effective magnetic field. Taking the 
derivatives of both sides of this equality with respect to scale "fc" yields, after rearrangement 

- d k U k (4>) + <t0kU' k (4>) + \d k e k <t> 2 = -\d k e k < ip 2 > + {d k U' k {<f>) + <f>d k e k ) < tp > (A.2) 



Note that we have 



< ip > — <f> — 9 In z/dh 
< ip 2 > - < <p > 2 = d4>/dh = l/(dh/d<f>) 

1 



Inserting these identities into (|A.2[) yields 



(A.3) 



2 + 

A further simplification occurs with the choice e k oc k 2 which entails 



dM4>) = 1 T JT, ■ (A.4) 
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TABLE I: Critical parameters of the <3? 4 theory on a 3D simple cubic lattice in the LPA approximation using an USCO cut-off. 
From left to right : ~g, r c (g), n c , and A c . The data were obtained by fixing g and determining r c (g) by dichotomy. k c and A c 
were then obtained from (g~, r c (g)) via Eqs. ([5]). An uncertainty of ±f affects the last figure. 



9 


r c 




A c 


9 


r c 




Ac 


0.000000 


0.0 


0.166667 


0.0 


0.200000 10^ 


-0.191553 10 1 


0.187467 


0.117147 


0.100000 


-0.125538 lO" 1 


0.166861 


0.464044 10 _a 


0.250000 10^ 


-0.232170 10 1 


0.190042 


0.150484 


0.200000 


-0.249786 10" 1 


0.167052 


0.930213 10" a 


0.300000 10^ 


-0.271343 10 1 


0.192047 


0.184411 


0.300000 


-0.372955 10" 1 


0.167240 


0.139846 10^ 


0.350000 10^ 


-0.309339 10 1 


0.193600 


0.218640 


0.400000 


-0.495161 10" 1 


0.167425 


0.186875 lO - * 


0.400000 10^ 


-0.346354 10 1 


0.194792 


0.252959 


0.500000 


-0.616484 10 _1 


0.167608 


0.234105 10~ a 


0.450000 10^ 


-0.382536 10 1 


0.195693 


0.287219 


0.600000 


-0.736986 lO" 1 


0.167789 


0.281532 10 _ * 


0.500000 10^ 


-0.418000 10 1 


0.196361 


0.321312 


0.700000 


-0.856719 lO^ 1 


0.167968 


0.329154 lb'' 2 - 


0.550000 10^ 


-0.452842 10 1 


0.196839 


0.355168 


0.800000 


-0.975725 10" 1 


0.168144 


0.376968 10~^ 


0.600000 10^ 


-0.487136 10 1 


0.197164 


0.388737 


0.900000 


-0.109404 


0.168319 


0.424970 VT* 


0.650000 10^ 


-0.520949 10 1 


0.197365 


0.421990 


0.100000 10 1 


-0.121170 


0.168492 


0.473160 V)-' 1 


0.700000 10^ 


-0.554300 10 1 


0.197458 


0.454881 


0.120000 10 1 


-0.144516 


0.168833 


0.570091 10~^ 


0.750000 10^ 


-0.587340 10 1 


0.197484 


0.487499 


0.140000 10 1 


-0.167631 


0.169167 


0.667742 VT* 


0.800000 10^ 


-0.620008 10 1 


0.197437 


0.519751 


0.160000 10 1 


-0.190530 


0.169495 


0.766098 10 -a 


0.850000 10^ 


-0.652375 10 1 


0.197337 


0.551677 


0.180000 10 1 


-0.213227 


0.169818 


0.865143 10~^ 


0.900000 10^ 


-0.684407 10 1 


0.197183 


0.583218 


0.200000 10 1 


-0.235735 


0.170135 


0.964863 VT* 


0.950000 10^ 


-0.716333 10 1 


0.197020 


0.614599 


0.300000 10 1 


-0.345763 


0.171648 


0.147315 10 _i 


0.100000 ltf 


-0.747982 10 1 


0.196819 


0.645628 


0.400000 10 1 


-0.452243 


0.173055 


0.199654 10 _i 


0.110000 10 J 


-0.810742 10 1 


0.196364 


0.706911 


0.500000 10 1 


-0.555797 


0.174373 


0.253384 lO^ 1 


0.120000 10^ 


-0.872935 10 1 


0.195870 


0.767298 


0.600000 10 1 


-0.656858 


0.175612 


0.308396 10" 1 


0.130000 W 


-0.934719 10 1 


0.195365 


0.826962 


0.700000 10 1 


-0.755749 


0.176780 


0.364598 10" 1 


0.140000 10 J 


-0.996241 10 1 


0.194871 


0.886080 


0.800000 10 1 


-0.852719 


0.177884 


0.421903 10 _i 


0.150000 10^ 


-0.105764 10^ 


0.194406 


0.944839 


0.900000 10 1 


-0.947970 


0.178929 


0.480232 10 _i 


0.160000 10 J 


-0.111905 10^ 


0.193981 


0.100343 10 1 


0.100000 10^ 


-0.104167 10 1 


0.179919 


0.539513 10" 1 


0.170000 10 J 


-0.118061 10^ 


0.193609 


0.106206 10 1 


0.150000 10^ 


-0.149095 10 1 


0.184166 


0.847930 lO" 1 


0.172979 10^ 


-0.119900 10^ 


0.193509 


0.107956 10 1 
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TABLE II: Critical parameters of the <3? 4 theory on a 3D simple cubic lattice in the LPA approximation using an LMD cut-off. 
From left to right : ~g, r c (g), n c , and A c . The data were obtained by fixing g and determining r c (g) by dichotomy. k c and A c 
were then obtained from (g~, r c (g)) via Eqs. ([5]). An uncertainty of ±f affects the last figure. 



9 


r c 




A c 


9 


r 




A c 


0.000000 


0.0 


0.166667 


0.0 


0.800000 10 2 


-0.625341 10 1 


0.198459 


0.525147 


0.100000 


-0.125603 10" 1 


0.166861 


0.464045 10 -3 


0.850000 10 2 


-0.656974 10 1 


0.198190 


0.556459 


0.200000 


-0.250055 10" 1 


0.167053 


0.930221 10 -3 


0.900000 10 2 


-0.688192 10 1 


0.197864 


0.587250 


0.300000 


-0.373541 10" 1 


0.167242 


0.139849 10^ 


0.950000 10 2 


-0.719022 10 1 


0.197488 


0.617527 


0.400000 


-0.496161 10 _1 


0.167428 


0.186881 10 -a 


0.100000 10 3 


-0.749487 10 1 


0.197074 


0.647300 


0.500000 


-0.617986 10" 1 


0.167613 


0.234117 10~ 2 


0.110000 10 3 


-0.809410 10 1 


0.196151 


0.705380 


0.600000 


-0.739072 10" 1 


0.167795 


0.281552 10~ 2 


0.120000 10 3 


-0.868108 10 1 


0.195140 


0.761593 


0.700000 


-0.859458 10" 1 


0.167976 


0.329184 10 -2 


0.130000 10 3 


-0.925704 10 1 


0.194072 


0.816050 


0.800000 


-0.979189 10 _1 


0.168154 


0.377011 10~ 2 


0.140000 10 3 


-0.982301 10 1 


0.192969 


0.868861 


0.900000 


-0.109830 


0.168331 


0.425031 10~ 2 


0.150000 10 3 


-0.103800 10 2 


0.191849 


0.920148 


0.100000 10 1 


-0.121680 


0.168507 


0.473241 10~ 2 


0.160000 10 3 


-0.109284 10 2 


0.190719 


0.969966 


0.200000 10 1 


-0.237346 


0.170181 


0.965382 HT* 


0.170000 10 3 


-0.114700 10 2 


0.189603 


0.101856 10 1 


0.250000 10 1 


-0.293536 


0.170971 


0.121796 10 _i 


0.172979 10 3 


-0.116300 10 2 


0.189272 


0.103280 10 1 


0.300000 10 1 


-0.348802 


0.171735 


0.147464 10" 1 


0.180000 10 3 


-0.120028 10 2 


0.188475 


0.106569 10 1 


0.400000 10 1 


-0.456919 


0.173190 


0.199965 10" 1 


0.190000 10 3 


-0.125300 10 2 


0.187373 


0.111177 10 1 


0.500000 10 1 


-0.562248 


0.174560 


0.253927 10" 1 


0.200000 10 3 


-0.130506 10 2 


0.186282 


0.115670 10 1 


0.600000 10 1 


-0.665180 


0.175854 


0.309247 10 _i 


0.225000 10 3 


-0.143280 10 2 


0.183644 


0.126469 10 1 


0.700000 10 1 


-0.766006 


0.177080 


0.365834 10" 1 


0.250000 10 3 


-0.155743 10 2 


0.181139 


0.136714 10 1 


0.800000 10 1 


-0.865000 


0.178243 


0.423608 lO" 1 


0.275000 10 3 


-0.167936 10 2 


0.178772 


0.146480 10 1 


0.900000 10 1 


-0.962205 


0.179345 


0.482472 10 _i 


0.300000 10 3 


-0.179892 10 2 


0.176537 


0.155827 10 1 


0.100000 10 2 


-0.105791 10 1 


0.180395 


0.542370 10" 1 


0.350000 10 3 


-0.203196 10 2 


0.172445 


0.173467 10 1 


0.150000 10 2 


-0.151715 10 1 


0.184929 


0.854965 10" 1 


0.400000 10 3 


-0.225821 10 2 


0.168797 


0.189951 10 1 


0.200000 10 2 


-0.195105 10 1 


0.188483 


0.118420 


0.450000 10 3 


-0.247888 10 2 


0.165533 


0.205508 10 1 


0.250000 10 2 


-0.236560 10 1 


0.191267 


0.152429 


0.500000 10 3 


-0.269479 10 2 


0.162590 


0.220296 10 1 


0.300000 10 2 


-0.276459 10 1 


0.193433 


0.187082 


0.550000 10 3 


-0.290700 10 2 


0.159940 


0.234491 10 1 


0.350000 10 2 


-0.315063 10 1 


0.195101 


0.222043 


0.600000 10 3 


-0.311562 10 2 


0.157522 


0.248133 10 1 


0.400000 10 2 


-0.352561 10 1 


0.196364 


0.257059 


0.650000 10 3 


-0.332136 10 2 


0.155317 


0.261337 10 1 


0.450000 10 2 


-0.389103 10 1 


0.197299 


0.291951 


0.700000 10 3 


-0.352456 10 2 


0.153295 


0.274161 10 1 


0.500000 10 2 


-0.424800 10 1 


0.197965 


0.326583 


0.750000 10 3 


-0.372553 10 2 


0.151435 


0.286658 10 1 


0.550000 10 2 


-0.459733 10 1 


0.198407 


0.360850 


0.800000 10 3 


-0.392455 10 2 


0.149718 


0.298872 10 1 


0.600000 10 2 


-0.494000 10 1 


0.198672 


0.394704 


0.850000 10 3 


-0.412184 10 2 


0.148127 


0.310840 10 1 


0.650000 10 2 


-0.527631 10 1 


0.198782 


0.428072 


0.900000 10 3 


-0.431761 10 2 


0.146650 


0.322593 10 1 


0.700000 10 2 


-0.560707 10 1 


0.198771 


0.460949 


0.950000 10 3 


-0.451203 10 2 


0.145275 


0.334160 10 1 


0.750000 10 2 


-0.593264 10 1 


0.198657 


0.493310 


0.100000 10 4 


-0.470524 10 2 


0.143992 


0.345561 10 1 
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TABLE III: Critical line k c (A) of the $ 4 model on a 3D simple cubic lattice. The MC data for k c (second column) are those of 
Ref. [10( 1 . The LPA data reported in the third column were obtained with the use of LMD regulator by fixing A and determining 
k c (A) by dichotomy. The (signed) relative error is defined as e = (k c ,lpa — k c ,mc)/k c ,mc. 



A 


Kc.MC 


K C ,LPA 


e 


0.1 


0.18670475 


0.1866196 


-0.000456 


0.2 


0.19421255 


0.1941031 


-0.000564 


0.4 


0.19879185 


0.1986966 


-0.000479 


0.7 


0.19626510 


0.1962421 


-0.000117 


0.8 


0.19438785 


0.1943946 


0.000003 


0.9 


0.19225565 


0.1922928 


0.000193 


1.145 


0.18644630 


0.1865659 


0.000641 


1.3 


0.18261165 


0.1827799 


0.000921 


1.4 


0.18013945 


0.1803395 


0.001111 


1.5 


0.17769270 


0.1779256 


0.001311 


2.5 


0.15671735 


0.1572028 


0.003098 
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